Magnetic Bose glass phases of coupled antiferromagnetic dimers with site dilution 
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We numerically investigate the phase diagram of two-dimensional site-diluted coupled dimer systems in an 
external magnetic field. We show that this phase diagram is characterized by the presence of an extended Bose 
glass, not accessible to mean-field approximation, and stemming from the localization of two distinct species 
of bosonic quasiparticles appearing in the ground state. On the one hand, non-magnetic impurities doped into 
the dimer-singlet phase of a weakly coupled dimer system are known to free up local magnetic moments. The 
deviations of these local moments from full polarization along the field can be mapped onto a gas of bosonic 
quasiparticles, which undergo condensation in zero and very weak magnetic fields, corresponding to transverse 
long-range antiferromagnetic order. An increasing magnetic field lowers the density of such quasiparticles to 
a critical value at which a quantum phase transition occurs, corresponding to the quasiparticle localization on 
clusters of local magnets (dimers, trimers, etc.) and to the onset of a Bose glass. Strong finite-size quantum 
fluctuations hinder further depletion of quasiparticles from such clusters, and thus lead to the appearance of 
pseudo-plateaus in the magnetization curve of the system. On the other hand, site dilution hinders the field- 
induced Bose-Einstein condensation of triplet quasiparticles on the intact dimers, and it introduces instead a 
Bose glass of triplets. A thorough numerical investigation of the phase diagram for a planar system of coupled 
dimers shows that the two above-mentioned Bose glass phases are continuously connected, and they overlap in 
a finite region of parameter space, thus featuring a two-species Bose glass. The quantum phase transition from 
Bose glass to magnetically ordered phases in two dimensions is marked by novel universal exponents (2 ~ 2, 
P « 0.9, v ~ 1). Hence we conclude that doped quantum antiferromagnets in a field represent an ideal setting 
for the study of fundamental dirty-boson physics. 

PACS numbers: 75.10.Jm, 72.15.Rii, 05.30.Rt, 03.75.Lm 



I. INTRODUCTION 

Quantum phase transitions in spin-gapped antiferromagnets 
(AF) have been extensively studied both theoretically and ex- 
perimentally. Several mechanisms exist which can drive these 
systems from a quantum disordered ground state with a finite 
spin gap to a gapless antifeiTomagnetic long-range ordered 
(AFLRO) state. Possible ways to close the spin gap are the 
modulation of the superexchange couplings by, e.g., mechani- 
cal deformation of the crystal under pressure', the application 
of a magnetic field2ri2ii^ii2 or the dilution of the system with 
non-magnetic impuritie s"i^"^''^ . 

For definiteness we will discuss the case of spin-gapped 
systems with a singlet ground state, associated with the mag- 
netic Hamiltonian of weakly coupled dimers. In these systems 
a finite singlet-triplet spin gap opens for sufficiently weak 
inter-dimer couplings. The application of a magnetic field 
to the system can close the gap at a critical value at which 
a quantum phase transition occurs, which is well described as 
a Bose-Einstein condensation (BEC) of 5 = 1 triplets^ii^— . 
The field has the combined effects of partially polarizing the 
dimers and inducing transverse AFLRO, which corresponds 
to a superfluid (SF) state of the field-induced triplets. Alterna- 
tively, doping weakly-coupled dimer systems with static non- 
magnetic impurities also has dramatic effects. Each dopant 
releases a spin- 1/2 degree of freedom out of the singlet, and 
this S = 1/2 moment gets exponentially localized around 
the site of the spin left unpaired^^. Hereafter this localized 
S = 1/2 degree of freedom will be denoted as local mo- 



ment (LM). The overlap of the localized wave functions of 
two nearby LMs produces effective couplings that decay ex- 
ponentially with the distance between the two impurities-'*, 
and which alternate in sign depending on whether the LMs 
belong to the same sublattice or not, and consequently they 
are not frustrated. The resulting random network of LMs 
therefore develops AFLRO at T = for any finite doping 
concentrations^,, and its low-energy excitations give rise to a 
phenomenon of order-by-disorder (OBD)^^'^^. 

It is intriguing to investigate the nature of the ground state 
in a spin-gapped system when both factors perturb the gapped 
state, i.e., a magnetic field and non-magnetic impurities, are 
present. In this paper we focus on the interplay between 
site dilution and an applied magnetic field in a planar cou- 
pled dimer system. We provide a detailed picture of the rich 
phase diagram of the system, with particular emphasis on the 
emergence of an extended novel quantum-disordered phase 
with the nature of a Bose glass (BG), as introduced for three- 
dimensional systems in Refs. 28,29 and for two-dimensional 
systems in Refs. i30l432l In particular we concentrate on the 
appearance of a disordered-local-moment (DLM) phase ap- 
pearing in the diluted system upon application of a weak mag- 
netic field^^'^^. This phase is characterized by a highly inho- 
mogeneous response of the network of LMs, with a portion of 
the network (made of nearly isolated LMs) being completely 
polarized by the field, and another substantial remaining por- 
tion of interacting LM clusters which takes a partially polar- 
ized state, namely it hosts localized magnetic "holes" in the 
background of LMs aligned with the field. Due to the dis- 
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Crete structure of the clusters supporting the holes (namely 
dimers, trimers etc. of LMs), the magnetization process in the 
DLM phase proceeds in smooth steps: its characteristic fea- 
ture is represented by pseudo-plateaus (PPs) associated with 
the magnetic response of small LM clusters subject to a broad 
distribution of local fields exerted on them by the remaining 
polarized LMs. 

We here list the main results of this paper. 

1) We quantitatively link the features of the magnetization 
curve to the statistics of the LM clusters and the effective local 
fields. We compare the results for a site-diluted coupled-dimer 
system with those for an effective model of a random network 
of LMs, both in the quantum case of 5* = 1/2 and in the clas- 
sical (S = oo) limit. The random network of quantum LMs 
characterizes the DLM phase, while the comparison with the 
classical model highlights the strong quantum nature. 

2) We introduce a spin-to-boson mapping which leads to 
a straightforward interpretation of the magnetic holes hosted 
by the LM clusters as localized bosonic quasi-particle states, 
whose population can be changed by infinitesimal changes of 
the applied field (namely the quasi-particle chemical poten- 
tial). Hence the resulting DLM phase is a compressible disor- 
dered phase, analogous to a Bose glass. 

3) We map out the entire phase diagram of the system for 
various inter-dimer couplings, showing strong deviations from 
the mean-field phase diagram of Ref. 34. The phase diagram 
shows that the DLM phase is indeed continuously connected 
to a higher field BG phase. A detailed quantum-critical scal- 
ing analysis reveals consistent critical exponents for the OBD- 
to-DLM transition and for the BG-to-SF transition. Hence we 
provide numerical evidence of a novel universality class asso- 
ciated with the BG transition in two dimensions. 

The paper is organized as follows. In Sec.|ll]we describe the 
models investigated and the numerical tools used in this work. 
In Sec. Hn] we briefly review the numerical results exhibiting 
the sequence of quantum phases and quantum phase transi- 
tions induced by the field in the 2D site-diluted coupled-dimer 
system. In Sec.|lV]we discuss the structure of the magnetic 
response and energetics in the DLM phase, and we quantita- 
tively explain the main features of the DLM phase. In Sec. |V] 
we discuss the phase diagram of the 2D site-diluted dimer sys- 
tem as a function of inter-dimer coupling strengths in the pres- 
ence of a magnetic field. 



II. MODELS AND NUMERICAL METHOD 

We consider a two-dimensional 5 = 1/2 Heisenberg an- 
tiferromagnetic model of weakly coupled dimers on a square 
lattice. ^^'^^ The Hamiltonian reads 

{ij) ^dimers (/m) Gintcr— dimers 

- hJ2e,St, (1) 

i 

where e takes the values or 1 with probabilities p and (1 — p), 
modeling the doping of the system with non-magnetic impu- 
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FIG. 1: 2D coupled-dimer model with site dilution. The ellipses 
indicate the formation of local singlets on the intact dimers, open 
circles indicate vacancies, and arrows represent local moments. 



rities. The pairs (ij) run on the strong dimer couplings J, and 
the pairs {Im) on the weaker inter-dimer couplings J' - see 
Fig.[T] In SecHniand |lV]we first focus on the weak-coupling 
limit, i.e., J' <C J. Later, in Sec.|V]we investigate the entire 
phase diagram for J' < J. 

In the clean limit (p — 0) and at zero magnetic field there 
is a critical value {J' / J)c — 0.523 separating the magneti- 
cally ordered phase for J'/ J > {J' / J)^. from a dimer-singlet 
phase for J'/ J < [J' /J)cr'^ Hence, for J' < J the ground 
state is in a quantum disordered phase with a finite spin gap 
Aq « J — zJ' + 0( J'^), where z is the coordination number 
for the dimers. This gap can be overcome by a magnetic field 
he = Ao, at which a quantum phase transition to an antifer- 
romagnetically ordered phase occurs. This phase transition is 
well described by the picture of Bose-Einstein condensation 
of triplet quasiparticles.''''^ 

In a randomly doped system (p > 0), the zero-field ground 
state deviates strongly from a dimer singlet. As described in 
the introduction, any arbitrarily small (but finite) site dilution 
with concentration p creates an interacting random network of 
LMs which orders at T = 0.^^ When looking at the physics of 
the system for energies much below the gap Aq, and conse- 
quently for fields h <C Aq, the behavior of the diluted system 
Eq. ([Til can be captured by an effective spin-1/2 Heisenberg 
model for the random network of LMs, whose Hamiltonian 
reads 

H = Y,jlf^S.-S,-hJ2si (2) 

i<j i 

Here the S operators represent S' = 1/2 LMs, exponentially 
localized around the sites of unpaired spins. Hence, on a 
L X L lattice the number of effective spins is approximately 
Ns = pL^ (ignoring the case in which a dimer is fully re- 
placed by the impurities, and which does not produce any 
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LM). Since the diluted sites are randomly located throughout 
the system in Eq. ([TJ, in the effective model Eq. (|2]i the LMs 
are also assumed to be randomly distributed on the square lat- 
tice. The effective interactions Jij^^ between LMs take 
the asymptotic form for large inter-moment distance^l^^i^i 

j(eff)=(_)|.-,|-i:^,-n./Co^ (3) 

where ,/i is a parameter determined by J and J' of the origi- 
nal diluted model, ^ '^o^ '■^^ correlation length of the 
gapped ground state in the clean limit, jfj^'' has ferromag- 
netic /antiferromagnetic character depending on whether the 
sites i and j belong to the same/different sublattices; hence 
the couplings are not frustrated, and they favor AFLRO of the 
LMs. a is a dimension-dependent exponent, taking values 
a = (ID), a = 1 (2D), and a ^ 3/2 (3D). In the following 
we will concentrate on the 2D case. When rij = 1 we need 

to recover the limit Jij^"* = Jie~^/'>° = J' for the direct in- 
teraction of two unpaired spins through a J' coupling. This 
leads to the identification Ji = J'e^^^°, which in turn implies 

= (_)l-.|-iile-(n,-i)/?o. (4) 
rij 

Actually the second-order perturbation calculation which de- 
livers the asymptotic form of the effective couplings, Eq. 
gives Jl ~ {zJ')'^/J; this means that the choice Ji ^ J', 
which correctly reproduces the short-range behavior of the ef- 
fective couplings, leads to an overestimation of such couplings 
in the long-distance limit. Such an overestimation is nonethe- 
less irrelevant, given the exponential suppression of Jij^'' at 
large inter-moment distances. 

In the following, we apply the Stochastic Series Expansion 
(SSE) quantum Monte Carlo (QMC)" to both the original and 
effective Hamiltonians, given by Eq. [T]and |2] in order to in- 
vestigate the evolution of the ground state of the system as a 
function of the applied magnetic field. For the effective model 
with Ns randomly distributed spins and non-frustrated long- 
range exchange interactions, the SSE approach has the advan- 
tage of requiring a computational effort 0{Ns In A^^) instead 
of the naive 0{N'^)?^ To efficiently reach the ground-state 
behavior via the finite-temperature SSE we make use of a (3- 
doubling approach^^ up to an inverse temperature /3J = 2^^. 
The full statistics of the disorder distribution is well repro- 
duced by averaging all results over « 300 disorder realiza- 
tions. 

In the presence of an applied magnetic field, spontaneous 
AFLRO is reflected by the appearance of a finite transverse 
staggered magnetization, estimated as 

= ^5^K7r)/L2, (5) 

where L is the linear system size and S'^(7r,7r) is the trans- 



verse staggered structure factor, defined as 

S^i7r,n) = ±^J2i-f''^{SfSJ + SfSy). (6) 

The summation runs over all pairs of occupied sites in the 
original site-diluted model, and it runs over all possible pairs 
of LMs in the effective spin model. We also estimate the spin 
stiffness, corresponding to a superfluid density upon spin-to- 
boson mapping, via the fluctuations of the winding numbers 
Wx,Y of the SSE worldlines along the two lattice directions 
X and yiS: 

Ps^ T^j{W^ + W^)- (7) 



III. MONTE CARLO RESULTS FOR WEAKLY COUPLED 
DIMERS 

In this section we present the QMC results for the full 
evolution of the ground state of site-diluted weakly coupled 
dimers upon varying an applied magnetic field. Here we 
focus on the case of weakly coupled dimers by choosing 
J'/ J = 1/4. In the clean case, this system is far away from 
the critical point sa 0.523 J in zero field, and it exhibits 
a disordered state with a correlation length ~ 1- When a 
magnetic field is applied, the clean system experiences a quan- 
tum phase transition to a field-induced antiferromagnetically 
ordered phase at h'f^ — Aq = 0.60(1) J. Doping with non- 
magnetic impurities produces LMs. Our choice of J' <C Aq 
guarantees that the energy scale for the interactions between 
LMs is well separated from the energy scale of interactions 
for the intact dimers. Hence the application of a weak field 
h ^ J to the system will exclusively probe the response of 
LMs, while intact dimers remain essentially frozen in local 
singlets. 

The succession of field-induced phases in the site-diluted 
system is fully captured by investigating the staggered mag- 
netization, rus (Eq. |5]l, and the uniform magnetization, rriu — 
'^i{Sl)/LF'- The former detects the presence of order in 
the ground state, while the latter probes the low-energy spec- 
trum both in the ordered and in the disordered phase, exhibit- 
ing the presence of a triplet gap via magnetization plateaus. 
The field dependence of these quantities at a dilution con- 
centration p = 1/8 is shown in Fig.|2] The following suc- 
cession of phases is observed in the field range < /i < 
J: order-by-disorder phase (OBD), disordered-local-moment 
phase (DLM), plateau phase (PL), Bose glass phase (BG), 
XY-ordered phase (XY). 

At zero field, a finite staggered magnetization develops in 
the system due to the effective interactions among LMs: an 
extrapolation of the QMC data to infinite size gives — 
0.032(3) (see Fig.|2jb)), indicating that the system is in the 
OBD phase. Further careful finite-size scaling analysis, dis- 
cussed in Sec. [Vl shows that the OBD phase also survives at 
smafl finite fields up to /i = /idlm = 0.007(1) J < J'. At 
this critical field the uniform magnetization takes the value 
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FIG. 2: (a): Field dependence of the uniform magnetization, the uni- 
form susceptibility, and the staggered magnetization, of weakly cou- 
pled dimers. The coupling ratio is J'/ J — 1/4, and the dilution 
concentration is p = 1/8. The orange diamonds show the staggered 
magnetization in the thermodynamic limit, extrapolated by finite- 
size scaling as shown in (b). OBD: order-by-disorder phase; DLM: 
disordered-local-moment phase; PL: plateau phase; BG: Bose glass 
phase; and XY: XY-ordered phase (magnetic condensate). 



system. The reason for this is the finite (albeit exponentially 
small) probability of finding an arbitrarily large clean region 
which is devoid of any impurity; the local gap of such a clean 
region can be arbitrarily close to Aq so that a field h > Aq 
can create a localized dimer triplet in that region. At vari- 
ance with the clean system, the field-induced localized quasi- 
particles do not condense in an extended state with finite su- 
perfluid response (spin stiffness), but they are instead local- 
ized by disorder and give rise to a quantum disordered Bose 
glass state. ^"'^^ In fact the fully polarized LMs act essentially 
as impenetrable barriers to dimer triplets, while intact dimers 
close to the non-magnetic impurities have higher local gaps 
to dimer triplets than in the clean regions. Hence site dilution 
creates an effective disorder potential which has the effect of 
Anderson-localizing the dimer triplet quasiparticles appearing 
in the ground state^i 

The resulting BG phase extends over the field range 
0.60(1) h/J ^ 0.67(1) in Fig.l^a), and, similarly to the 
DLM phase, it is identified by the absence of a gap, as shown 
by the finite slope in the magnetization curve (namely a finite 
susceptibility Xu = dmu/ dh)^— 

At Hbg/J ~ 0.67, the applied magnetic field finally 
drives the system to an XY-ordered antiferromagnetic state. 
In the picture of bosonic triplet quasiparticles, this magneti- 
cally ordered state is a condensate with finite superfluid re- 
sponse. Hence this system realizes a magnetic BG-SF tran- 
sition, with critical exponents which are quite different from 
the more conventional Mott insulator to superfluid (MI-SF) 
transition. '"'"^^ As h is increased to reach half filling of the 
dimer triplets, the behavior of the system is better understood 
in terms of singlet quasiholes (given that the triplet quasiparti- 
cles are effectively hardcore). Such quasiholes experience the 
reverse succession of phases to that of triplet quasiparticles, 
namely a transition from XY to BG and then from BG to a 
plateau phase corresponding to full polarization of the entire 
sample.-''^ 



TO„ = 0.0208(6), namely it is much less than the value 
— = 1/16 corresponding to full saturation of the 
local moments. In fact, the saturation value is only attained at 
a much higher field h = /ipL ~ 2 J'. Hence the ground state 
of the diluted system at /idlm < h < /ipL has a striking fea- 
ture: the LMs have lost spontaneous ordering transverse to the 
field, but they are far from being all polarized. In contrast, in 
a clean antiferromagnet the destruction of a spontaneously or- 
dered state is always accompanied by the full polarization of 
its constituents. Moreover for h > /idlm the magnetization 
keeps growing with the field, revealing a gapless spectrum in 
the absence of spontaneous ordering. 

For /ipL < h < Aq the uniform magnetization remains 
constant at m^^*, indicating a gapped, disordered plateau 
phase (PL)^"* in which the field completely aligns all the un- 
paired spins, but is not sufficiently strong to break the singlets 
on strongly coupled dimers. The existence of the PL phase is 
then a characteristic of the weakly coupled dimer system, and 
we will see that it disappears when J' ^ Aq. 

When h ^ Aq the gap to excite a triplet in a region of in- 
tact dimers (dimer triplet) closes just as it does in the clean 



IV. FIELD-INDUCED DLM PHASE 
A. Magnetization curve and effective model 

In this section we focus on the unconventional properties 
of the DLM phase. As already pointed out in the previous 
section, this phase is characterized by the absence of a gap, 
as revealed by a finite susceptibility. This is actually not at 
all obvious from Fig. |2] where the susceptibility appears to 
nearly vanish in correspondence with intermediate pseudo- 
plateaus (PP) at 73% and 95% of the saturation magnetiza- 
tion. Only a careful analysis of the temperature scaling of the 
susceptibility^^ shows that this quantity remains finite at all 
field values in the DLM phase down to T = 0. The first PP 
extends up to ss 0.7J'; the second appears at /i w 1.2J'; 
and the true saturation plateau is only attained at /i w 2 J'. 

As argued in Sec. |ll]the physics of the site-diluted dimer 
model at low energy (namely at low fields for T = 0) can be 
fully captured by an effective model of coupled LM, Eq. 
Here we make this statement more quantitative, by directly 
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FIG. 3: Field dependence of ms, rriu, and Xu in the effective model, 
Eq. ([2}, with = 1.0 at a dilution concentration p — 1/8. 



investigating the physics of the effective model via QMC sim- 
ulations. The field evolution of the uniform and staggered 
magnetization for the model of Eqs. (|2|l and ([3]) with = 1 
is shown in Fig. [3] Similarly to the results for the original di- 
luted model, we observe the destruction of the AFLRO phase 
of the LMs at very small fields, and the occurrence of an ex- 
tended DLM phase. The occurrence of PPs in the magnetiza- 
tion curve is fully captured in the effective model, albeit with 
less sharp features. The full qualitative correspondence be- 
tween the original Hamiltonian and the effective LM model 
corroborates the picture that the magnetic response of the sys- 
tem up to the plateau phase is dominated by LMs. In the 
following subsection we will bring this observation to a mi- 
croscopic level, relating the occurrence of PPs to the quantum 
magnetic response of LM clusters. 



B. Quantum nature of the disordered-local-moment phase 

In order to emphasize the genuinely quantum character of 
the DLM phase, in this subsection we focus on a comparison 
of the system of quantum 5=1/2 LMs described by Eq. (|2]i 
with its classical limit 5 = oo. The classical limit of Eq. (|2]l 
is investigated via a classical Monte Carlo study of a 40 x 40 
lattice with a concentration p = 1/16 of LMs. The ground 
state of the classical system is reached by careful thermal an- 
nealing with a linear protocol, and results are averaged over 
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FIG. 4: (a): Classical S oo local moments form a canted antifer- 
romagnetic order at finite field, (b): In the quantum case, neighbor- 
ing local moments prefer to form strongly correlated states on local 
clusters, leading to the novel DLM phase. The black and red arrows 
distinguish spins located on different sublattices. 
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FIG. 5: Field dependence of mu and of the transverse spin com- 
ponent 5'^'^' in the classical effective spin model with = 1-0, 
p = 0.0625 and L = 40. The corresponding m„ and nis of the 
quantum model with the same model parameters are also shown for 
comparison. The staggered magnetization ms of the quantum model 
has been magnified by a factor of 100. 



typically 500 — 1000 configurations. 

In the classical limit at T = 0, the system of randomly dis- 
tributed LMs in a sufficiently weak magnetic field exhibits a 
canted antiferromagnetic ground state (Fig. |4ja)). Upon in- 
creasing the field, the most isolated LMs minimize their en- 
ergy by fully aligning with the field, but the remaining LMs 
preserve long-range order transverse to the field because of 
the direct long-range effective interactions, as clearly exhib- 
ited by Fig. |5] Long-range order is destroyed only when the 
field exceeds the polarization value of the largest connected 
LM cluster that can form in the system, thereby polarizing all 
the LMs. This means that, in the classical limit, the destruc- 
tion of AFLRO is always accompanied by the full saturation 
of the magnetization. 
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Fig. |5] compares the classical behavior with the one of the 
quantum system with the same concentration of LMs. The 
uniform magnetization curve shows the emergence of a PP 
also at the classical level, precisely at the same value as in 
the quantum case. This result is consistent with the picture 
(both classical and quantum) of PPs originating from strongly 
coupled LMs which resist to being fully polarized by the field. 
This picture relies on the strong inhomogeneity of the network 
of interacting LMs. This network exhibits an average distance 
between LMs of rave = 1/ \/v^ a corresponding average 
coupling Javo ^ J' in the limit of diluted LMs. Yet the mag- 
netic response is strongly dominated by the long tails of distri- 
bution of LM couplings and not by its average. In particular, 
there exist small clusters of LMs with linear size r <C r-ave, 
which interact with much stronger couplings j'^'^^^ ^ J' . 
The main examples of such clusters are nearest-neighboring 
dimers and trimers (see the cartoon of Fig. |4ja)). Evidently 
only a field h ^ J' can overcome the antiferromagnetic cor- 
relations developing in these clusters. 

In the classical limit, LM clusters not only retain short- 
range antiferromagnetic correlations due to their strong inter- 
nal coupling, but also long-range correlations due to the direct 
long-range couplings existing among them. Hence, AFLRO 
classically survives the strong field h ^ J', as shown in Fig.|5] 
In particular the transverse staggered magnetization shows a 
PP corresponding to that of the uniform magnetization, a fact 
which corroborates the above picture. 

On the contrary, the quantum mechanical behavior of the 
LM network is significantly different, as already shown by 
the data on transverse staggered magnetization in Fig. |2] At 
the quantum level the presence of strong short-range antiferro- 
magnetic correlations persisting on LM clusters does not co- 
exist with long-range antiferromagnetic correlations between 
different clusters. In fact a phenomenon of quantum clus- 
tering of antiferromagnetic correlations (dimerization, trimer- 
ization etc.) turns out to minimize the energy at the quan- 
tum level, as sketched in Fig. lUb). As it will be quantita- 
tively shown in the next subsection, LM dimers have the ten- 
dency to freeze in a singlet, LM trimers tend to freeze in a 
S'tot = 1/2 state aligned with the field, etc. LM dimer sin- 
glets are obviously uncorrected with one another, while LM 
trimers are also antiferromagnetically uncorrected because 
the Stot = 1/2 state does not reflect the same sublattice struc- 
ture for all trimers. Hence, except for a small field range 
slightly above zero field, the entire phase with partially po- 
larized LMs is quantum disordered. 



C. Local Moment Clusters 

The dominant role played by the strongly interacting clus- 
ters of LMs in the DLM phase clearly emerges from the distri- 
bution of the local magnetization {Sf), as well as of the pair 
correlations C{rij) = {Si ■ Sj), developed by the inhomoge- 
neous network of LMs. The quantum behavior of these small 
LM clusters is revealed again by comparison with the classical 
limit. 

The distribution of (5f ) in the classical model is shown in 




FIG. 6: Distribution of the z component of local moments Si in the 
effective spin model for classical spins with = 1-0, p = 0.0625, 
and L = 40. 



Fig. |6] Here the classical spin length is normalized to 1/2 to 
compare with the quantum case. At zero field, (5f ) is ob- 
tained by averaging the ordered ground state over all possible 
orientations. The corresponding probability distribution is: 



1 



(8) 



The presence of a very weak field further enhances two peaks 
located at (S'f ) = ±0.5. The peak at (5f ) = 0.5 corresponds 
to spins fully polarized by the magnetic field, whereas the 
peak at (S^) = —0.5 corresponds to spins antiferromagneti- 
cally coupled to the fully polarized ones (not shown in Fig.|6]l. 
For h/ J' > 0.1, an extra peak appears at a positive (S^), and 
it moves towards (S'f ) = 0.5 with increasing magnetic field: 
this peak corresponds to the partially polarized (canted) LMs. 
At the same time, the peak at (5f ) = —0.5 drops. This is 
fully consistent with the picture that classical spins are con- 
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tinuously polarized. 




FIG. 7: (a)-(e): Distribution of the z component of the local moments 
Si in the effective spin model for quantum spins in various fields 
with 5o = 1.0, p = 0.0625, and L = 40; (f): Field dependence of 
the uniform magnetization m„ for the same model parameters. 



The local magnetization histograms for the quantum system 
of LMs are shown in Fig. |7] and they exhibit significant dif- 
ferences with respect to the classical case, related to the quan- 
tum behavior of LM clusters. For ft, = (not shown in Fig.|7ll 
the histogram for (S'f ) exhibits a unique (5-peak at (S'f ) = 0, 
which obviously corresponds to the fact that the ground state 
of a quantum Heisenberg antiferromagnet on a finite-size sys- 
tem is a total singlet. The evolution of the local-magnetization 
histograms is strongly non-monotonic upon applying an in- 
creasing field, and it allows to precisely track the succession 
of phases in the LM network. In fact for a small field, corre- 
sponding to the OBD phase, the (S'f ) — peak evolves into 
a broad distribution of local magnetizations, which gives a di- 
rect insight into the microscopic structure of the OBD state. A 
significant portion of the LMs is in a canted antiferromagnetic 
phase with variable canting, depending on the local structure 



of LM couplings. This corresponds to the tail of the distribu- 
tion at (SI) > 0. At the same time a significant portion of 
the LMs is fully polarized by the field, as shown by the emer- 
gence of a peak at (Sf ) = 1 /2. These polarized LMs create an 
effective static field for the LMs which are coupled to them. 
In the case of antiferromagnetic couplings, this field points 
opposite to the external field, and it induces a fraction of the 
LMs to develop a negative local magnetization, as shown by 
the fat tail of the histogram at (5f ) < 0. Yet the persistence 
of a central peak at (5f) = shows that a large fraction 
of the LMs remains frozen in local singlets, associated with 
even-number LM clusters (dimers, quadrumers, etc.). At the 
same time two further peaks clearly emerge in the histogram, 
at (Sf ) = l/3andat {Sf) = —1/6: they can be easily associ- 
ated with the state |5triinor = 1/2, 5^ = 1/2) forLM trimers, 
where one has a trimer magnetization S"^ = X]ietrimor('^i^) — 
1/3 + 1/3 — 1/6 = 1/2. Further (and smaller) peaks are to 
be attributed to larger (and more rare) LM clusters. 

The singlet peak at (5f ) = persists over a significant field 
range. In particular in the field interval 0.3 < h/ J' < 0.8 
the local-magnetization histogram changes very little, corre- 
sponding to the appearance of the first PPs (as presented in 
Fig- Elf))- Only for h > J', namely when the field over- 
comes the largest possible triplet gap associated with a LM 
dimer-singlet, we observe a decrease of the singlet peak and 
correspondingly an increase in the slope of the magnetiza- 
tion curve, which allows then to associate the first PP pri- 
marily with the quantum magnetic response of LM dimers. 
Remarkably the LM-trimer peaks at (5f) — 1/3 and —1/6 
persist due to the larger gap to full saturation exhibited by 
these structure. The weak field dependence of such peaks for 
fields in the interval 1.2 < h/J' < 1.4 allows then to unam- 
biguously associate the second PP to the LM trimer gap from 

l^trimcr - 1/2, = 1/2) tO |^trimcr = 3/2, = 3/2), 

corrected by the local field interacting with each trimer in the 
network. 

The correspondence of the PP to the magnetic response of 
LM clusters is further confirmed by investigating the distri- 
bution of bond energies. The minimum bond energy is the 
one of the singlet Ei, — —jJ'- We hence define a normal- 
ized bond energy Et — — 4i?;,/(3 J'), upper-bounded by 1, 
and we plot the field dependence of the distribution of ln(£'b) 
in Fig. [8] In zero field we observe a very broad distribution, 
characteristic of a highly inhomogeneous system. In partic- 
ular there is a very broad peak at low energies (in absolute 
value), coiTesponding to the long tail of the distribution of LM 
couplings. But the most striking feature is a gap between the 
low-energy peak and two peaks corresponding to the stronger 
bonds, (Ef,) = 1 and (Ei,) =2/3. These peaks clearly cor- 
respond to dimer singlets and to trimer doublets, respectively. 
They are very resistant to the application of a field, while the 
low-energy part of the distribution changes continuously, de- 
veloping a strong peak which corresponds to the LMs fully 
polarized by the field. The dimer peak starts to decrease con- 
tinuously for h > 0.4J', and it disappears for h > J, cor- 
responding to the end of the first PP. This in turn reveals that 
the triplet gaps of the LM dimers follow a distribution, upper 
bounded by J', and due to the effect of the local effective field 
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FIG. 8: (a)-(d): Distribution of the normalized bond energies (Et) in 
the effective spin model with = 1-0, p = 0.0625 and system size 
L = 40. (e)-(h): Zoom on the distribution in the region exhibiting 
the dimer and trimer peaks. 



exerted by the other LMs surrounding the dimer and randomly 
coupled to it. Therefore the ensemble of LM dimers responds 
to a field in a continuous fashion, which gives rise to the fi- 
nite slope of the first PP. Similarly, the trimer peak starts to 
decrease only for h > 1.2, consistent with the appearance of 
the second PP with finite slope. 



tum Monte Carlo showing the persistence of Neel order. 

The effective model of interacting LMs, Eq. dU, shares sim- 
ilar features to those of a random-bond Heisenberg model, 
namely the presence of a broad distribution of energy scales 
for the LM couplings. As shown in this study, and as dis- 
cussed extensively in Ref.HH the system of LMs in zero field 
orders antiferromagnetically. We probe directly the possible 
formation of longer-range singlets by investigating the dis- 



tribution of the congelation function C{r 



- (S, ■ SA for 



antiferromagnetically coupled LMs at various inter-moment 
distances r^ . If sites i and j form an approximate sin- 
glet, then C{rij) w —3/4, while if they participate in a 
long-range trimer, C{rij) w —1/2. The results for rij = 
1, Vd, 3, VTS, Vl7 in a system with p = 1/16 are shown in 
Fig. |9] For r^j — 1, i.e. for nearest-neighboring moments, 
there are two sharp peaks in P{C{rij)) around C{rij) — 
—3/4 and C{rij) ~ —1/2 respectively. As discussed above, 
these two peaks correspond to LMs which participate in 
nearest-neighboring dimers and trimers. A similar peak struc- 
ture is also observed for — \/5 and rij — 3, indicating the 
appearance of similar structures at a larger length scale. The 
only difference is that the dimer and trimer peaks are weaker 
and broader, and that there is an extra peak around C(ry ) — 
for Tij = \/b and = 3, which corresponds to pairs of spins 
belonging to uncorrelated LM clusters {e.g. two short-range 
LM dimers or LM trimers). For > 3 the distribution func- 
tion completely changes its shape. There is a very broad peak 
around C{rij) ~ and a shoulder around C{rij) w —1/3, 
corresponding to long-range AF correlations, but no peaks 
at —3/4 or —1/2. It is interesting to notice that the aver- 
age inter- moment distance rave = 1/a/P = 4 roughly marks 
the boundary of two behaviors. Hence the above numerical 
findings imply that at zero field, LMs with inter-moment dis- 
tance rij < Tavo = 1/VP are generally involved in local sin- 
glets (doublets) associated with even (odd) clusters local, and 
they only marginally take part in the AFLRO. On the con- 
trary the long-range AF correlations are carried by LMs with 
rij }^ ?'avc- We finally notice that J^'^^^rij — ~ huLM- 
Thus, in the DLM phase almost all LMs involved in long- 
range clusters are fully polarized by the field, and thus the 
magnetic properties of DLM phase mostly depend on the 
statistics of the clusters made of nearest-neighboring LMs. 



D. Long-range clusters 

A fundamental result in one-dimensional quantum antifer- 
romagnets with random nearest-neighbor couplings is the for- 
mation of spin singlets at all energy scales and all distances, 
to form the so-called random-singlet phase. This phase is 
fully captured by a real-space renormalization group approach 
which, at each step, identifies the locally strongest antifeiTO- 
magnetic bonds and freezes the coiTesponding spins into a sin- 
glet, obtaining then effective couplings for the leftover spins 
via second-order perturbation theory around the singlet state. 
The application of a similar approach to the two-dimensional 
case of the Heisenberg antiferromagnet with random bonds 
fails to find a random-singlet phase,"*^ in agreement with quan- 



E. Statistics of Local Clusters and Pseudo-Plateaus: 
a Minimal Model 

As seen in the analysis of the QMC data presented in the 
previous section, the physics of the DLM phase can be mainly 
ascribed to the behavior of clusters of nearest-neighboring 
LMs immersed in the local field which is created by the sur- 
rounding isolated LMs. In this section we show how to cap- 
ture some fundamental aspects of the DLM phase with an el- 
ementary model inspired by the above results. 

The fundamental assumption we make is the diluteness of 
the LMs, namely p <C 1. Under this assumption the majority 
of LMs do not have nearest neighbors, and are typically at a 
distance 1 / ^ from the other LMs. We will denote these LMs 
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clusters are fully polarized, {S 
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C(rJ=<S^S>) 

FIG. 9: Distribution of inter-moment correlation functions C{rij) 
for various inter-moment distances Vij in the effective model with 
^0 = 1-0, p = 0.0625, and at zero field. Here sites i and j are 
always located on opposite sublattices. 



as isolated LMs. Assuming to have an applied magnetic field 
which is able to polarize most of the isolated LMs, we can es- 
sentially neglect the fluctuations of their magnetizations, and 
treat them as static variables. The remaining minority of LMs 
is arranged in clusters (dimers, trimers, etc.) which are rare 
and hence they seldom happen to be close to another. It is 
therefore a reasonable approximation to neglect all interac- 
tions between such clusters. Given that the LMs belonging to 
the clusters are the only dynamical variables left after freezing 
the isolated LMs, we can then model the magnetic response 
of the LM network in the DLM phase as the response of inde- 
pendent LM clusters immersed in the joint external field and 
in the local field created by the surrounding isolated LMs. For 
each cluster C we therefore need to diagonalize the following 
effective Hamiltonian 



c 



(cff) 



St (9) 



.^Qi = 1/2. This is of course 
not exact due to the existence of other LM clusters, but if they 
are sufficiently far apart from the cluster under consideration 

(which is true if p <C 1) we can neglect their couplings j'lf^'' 

to the cluster C, given the exponential decay of J^^^^ with the 
distance. 

After diagonalizing T-Lc for each cluster, we may calculate 
the uniform magnetization via 



{C} 



(10) 



where C represents the complement to the cluster C. In this 
simplified model we have assumed that all spins outside the 



where P{C) is the probability cluster C appears. Since for 
each finite cluster m„ (C) is a multi-step function at low tem- 
perature, m,j is expected to exhibit a multi-step character 
modulated by the distribution P{C). 

In Fig. [TOl' a) the magnetization curve for the effective 
model with = 1.0, p = 1/8, and L = 100 is calculated 
by applying Eq. [TO] We see that reaches the saturated 
value m^"* only when h/J' > 2. For h/J' ^ 2, the curve 
resembles the one obtained by QMC, i.e., several PPs can 
be resolved. By comparing with the magnetization curves of 
small clusters presented in Fig.[TOl"b)-(e), we can identify the 
first PP at h/J' < 0.7 to be associated with the 5* = sin- 
glet states of even-number clusters and the S" = 1/2 states 
of odd-number clusters. Since P{C) decays exponentially 
with the cluster size, this PP mostly comes from the S ~ 
state of LM dimers. Similarly we identify that the PP from 
1.0 < h/J' < 1.5 is associated with the 5=1/2 state of LM 
trimers. These are consistent with results obtained in previ- 
ous sections. Interestingly, we can resolve other weak struc- 
tures in the magnetic curve. For example, the little kinks at 
h/J' « 0.7 and h/J' « 1.7 are associated with the 5 = to 
S — 1 and S = Ito S — 2 transitions of four-site chains. We 
also notice that m„ saturates for h/J' ^ 2 for small clusters 
at least with sizes up to iV = 4. Actually, iV > 4 LM clusters 
are extremely rare so that they have only negligible contribu- 
tion to the magnetization curve. Hence we may take h — 2J' 
as the saturation field. 

We then compare the magnetization curve obtained by di- 
agonalizing LM clusters and then applying Eq. [TO] with the 
one via QMC for the same effective model. The results are 
presented in Fig.[TT]for two values of dilution concentration. 
We see that the curves obtained by the two different methods 
agree quite well, especially for the heights and positions of 
PPs. This agreement confirms that the PPs are ascribed to the 
distribution of small LM clusters. We notice that the curves 
from the two methods do not match well for very small h, 
and at h/J' ss 1. These discrepancies are closely related to 
the two assumptions for the LM clusters: since we assume all 
the single-site LMs are fully polarized by the magnetic field, 
the first PP calculated by the diagonalization of clusters is ex- 
pected to appear at very low field. Moreover the interactions 
among LM clusters, which were neglected in the diagonaliza- 
tion method, lead to smoother transitions between PPs, and to 
a significant smearing of the second PP. From Fig. [TTT b). we 
see that these features are captured by the QMC simulations. 
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FIG. 11: Field dependence of the normalized magnetization 
rnu/rn^^ in the effective model with ^ = 1.0 and p — 1/16 (in 
(a)) and p — 1/8 (in (b)). The results obtained by Eo.llOKred sym- 
bols) for a L — 100 lattice are compared with those from QMC 
simulations (black symbols) for a 1/ = 40 lattice. 



FIG. 10: Field dependence of the normalized magnetization 
m„/m^"* calculated by applying Eq.[T0]in the effective model with 
^0 = 1-0, p = 1/8, and L — 100 shown in (a). The normalized 
magnetization on small clusters up to A'^ = 4 with the same model 
parameters are also shown in (b)-(e). 



but they are missing in the diagonalization method. 



with parameters 



(off). 



= J, 



(cfT) 



1 



(12) 
(13) 
(14) 



F. Bose glass nature of the DLM phase 

As seen in the previous subsection, some of the most rel- 
evant features of the magnetization curve in the DLM phase 
can be captured by a model of independent clusters which are 
storing some LMs not aligned with the field. Upon a standard 
spin-boson transformation for S — 1/2 spins, which maps 
t/own-spins onto bosons (hereafter called LM quasiparticles) 
and M/?-spins onto bosonic holes (5*+ 6, 5^, and 

S"^ 1/2 - b^fb), the LM Hamiltonian, Eq. Q, takes the 
following form 



H^-J2u,iblb, 



K'j 



(11) 



The picture of the DLM phase as dominated by strongly 
quantum fluctuating LM clusters can be then recast in the 
bosonic language as a phase in which LM quasiparticles are 
expelled almost completely from the isolated sites, and they 
remain localized on the clusters, namely they resonate be- 
tween two sites on a LM dimer, they localize around the cen- 
tral site of a LM trimer, etc. We then apply the cluster decom- 
position of the Hamiltonian introduced in the previous sub- 
section, and we neglect inter-cluster couplings in the dilute 
limit. Neglecting the coupling between clusters and isolated 
spins as well, we obtain a ground state with isolated sites com- 
pletely empty of bosons, and clusters hosting one or more bo- 
son. Transfering a bosons from a cluster C to an isolated site 
has a chemical potential price Afj, ~ Ac ^ J' where is 
the gap to full polarization (or boson depletion) of the clus- 
ter. Given that all couplings between the cluster and the iso- 
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lated spins are smaller than J' we can introduce them pertur- 
batively. In doing so, we can still neglect the density-density 
coupling VijUiUj with i £ C and j e C, because the bosonic 
occupation of the isolated spins is negligible. The hopping 
terms ^ tij for i £ C and j G C lead in turn to a pertur- 
bative correction to the wave function of LM quasiparticles, 
which acquires an exponentially decaying tail (with localiza- 
tion length ^ I ln( Javc/ >^')l ^^) o^^r the nearest isolated sites, 
as it usually happens for a particle in a potential well of depth 
-Ac. 

The resulting picture of exponentially localized LM quasi- 
particles around the LM cluster locations is hence that of 
a Base glass. The density of states of the exponentially 
localized bosonic states is continuous, and correspondingly 
an infinitesimal change of the applied field always leads to 
adding/eliminating one boson, which is a characteristic fea- 
ture of a Bose glass. This picture will become evident when 
discussing the complete phase diagram of the system, as 
shown in the following section. 



V. PHASE TRANSITIONS IN THE 2D SITE-DILUTED 
COUPLED DIMER SYSTEM 

In this section, we discuss the complete T — phase di- 
agram of dimer systems with general coupling ratios ^ 
J'/ J ^ 1. This will allow us to give a unified picture of the 
novel quantum-disordered phases introduced by site dilution 
in the weakly coupled dimer systems. The main result is that 
the DLM phase, discussed in details in the previous section, 
is continuously connected to the Bose glass phase at higher 
fields. This is further evidence of the Bose glass nature of the 
DLM phase^^ 



A. Phase Diagram with Site Dilution 

We start by briefly recalling the T = phase diagram 
of a coupled-dimer system in the clean limit. At zero field, 
the ground state is quantum disordered and gapped for weak 
inter-dimer couplings, and it turns into a Neel-ordered state 
when the ratio J'/ J crosses the quantum critical point at 
( J'/ J)c ~ 0.523.-'^ This critical point belongs to the uni- 
versality class of the classical three-dimensional Heisenberg 
model^^'^s ^j^j^ ^ ^ 0.71(1), /3 = 0.36(1), and z = 1. In the 
presence of a magnetic field, the critical point at zero field ex- 
tends to a critical curve in the 2D [h/ ,])-{,]'/ J) plane (see 
dashed line in Fig. [T2] |. This curve separates the quantum 
disordered phase at low field and low inter-dimer coupling 
from the AFLRO phase. For this two-dimensional system, 
the universality class of this transition is the mean-field one 
(i^ = 1/2, ^ = 1/2, etc.) with z = 2.^^^ 

Fig. [12] shows the complete phase diagram of the two- 
dimensional doped system with p ^ 1/8. As already men- 
tioned in the introduction, doping with a small concentra- 
tion of non-magnetic impurities immediately induces antifer- 
romagnetic order in the zero-field ground state for every finite 
inter-dimer coupling J'. This is in agreement with the mean- 
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FIG. 12: Phase diagram of a site-diluted coupled dimer system at a 
dilution concentration p — 0.125. Three distinct phases are shown: 
the gapped plateau (PL) phase, corresponding to a bosonic Mott insu- 
lator; the disordered-local-moment (DLM) phase, corresponding to 
a Bose glass of LM quasiparticles; Bose glass (BG) phase of dimer 
triplets, and the two-BG phase, corresponding to a coexistence of 
LM-quasiparticle and dimer-triplet BG; and the XY ordered phase, 
corresponding to a condensate of magnetic quasiparticles. The phase 
boundary in the clean limit p — (black dashed line) is also shown 
for comparison. 



field approach of Ref. 34. For a large interval of J'/ J < 0.45 
(comprising the case J' = J/4 of Sec. Ullb . the system is 
driven back to a quantum disordered phase upon application 
of a field h = ft-oLM ^ J'. At a higher field /ibg (> Aq), 
the system experiences a second quantum phase transition 
back to the ordered phase, as discussed in Refs. [30l - [32i For 
J'/J < 0.275 the disordered phase is actually a multiple 
one, composed of two gapless regions - a DLM phase at low 
fields and a BG phase at higher field^— - and an interme- 
diate gapped plateau phase. The plateau phase is bounded 
from above by the critical line of the clean system, as expected 
from the fact that below that line no dimer triplet can appear 
on the rare regions of intact dimers. The curve that bounds 
the plateau phase from below, on the contrary, corresponds 
to ~ 2,/'; this is also to be expected as the lower critical 
field for this phase corresponds to the polarization field for all 
LM clusters whose inter-LM coupling is J', and whose ef- 
fective polarization field is effectively given by 2 J' (see also 
the discussion in Sec. IIV El and in Ref. 47). Hence the two 
curves bounding the plateau phase from above and from be- 
low have opposite slopes, and therefore cross each other at 
J'/J < 0.275, where the plateau phase disappears from the 
phase diagram, leaving space to a unique, continuous gapless 
disordered phase for 0.275 < J'/J < 0.45. 

Finally, for J'/J > 0.45 the low-field OBD phase merges 
with the high-field superfluid phase, marking the disappear- 
ance of all quantum-disordered phases in the system. It is 
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remarkable to observe that this happens well below the crit- 
ical value J'/J — 0.523 which marks the disappearance 
of the quantum-disordered phase in the clean limit. Hence 
site dilution of the lattice has the effect of significantly shift- 
ing the critical J'/ J ratio for the existence of a disordered 
phase in the system, as observed experimentally in Mg-doped 
TlCuCl34^ 

The topology of the phase diagram of Fig. [12] is generic 
for lattices of weakly coupled antiferromagnetic dimers (with- 
out frustration) with a concentration of vacancies that falls 
below the percolation threshold of the lattice. In particular 
this phase diagram will apply both XalD and to 31? arrays 
of dimers, which are percolating up to a finite concentration 
of vacancies. In particular a fundamental feature is that the 
disordered and gapless phases (DLM, 2BG and BG) are sim- 
ply connected in the phase diagram, so that they systemati- 
cally separate the gapped phase of fully polarized LMs from 
the magnetically ordered phases. In other words, in presence 
of site dilution there is no direct transition between a spin- 
gapped phase and a magnetically ordered phase. This par- 
ticular aspect of the topology of the phase diagram is a gen- 
eral property of that of bosons in a disordered lattice'*'-'*'', with 
the correspondence between spin-gap phases and Mott-/band- 
insulating phases, and between magnetically ordered phases 
and superfluid phases. 

When considering vacancy concentrations larger than that 
studied here (but still below the percolation threshold), one 
expects the gapless phases (DLM and Rose glass) to grow at 
the expenses of the plateau phase. The magnetically ordered 
phase at low fields (OBD phase) will also grow towards larger 
fields, as a denser network of LMs develops a magnetic or- 
der which is more robust to the application of a field. On the 
other hand, the ordered phase at large fields will instead shrink 
at the expenses of the Bose glass phase, given that a larger 
triplet density (namely a larger magnetization) is necessary 
for the triplet gas to condense when the disorder is higher It 
is important to contrast the phase diagram of Fig. [12] with the 
one presented in Ref. ,34 on the basis of mean-field calcula- 
tions. In the mean-field phase diagram all the gapless disor- 
dered phases, dominated by Anderson localization of triplets 
and/or moments aligned oppositely to the field, are completely 
absent. This is easily understandable since Anderson local- 
ization cannot be captured at the mean-field level; the gapless 
nature of the above mentioned phases comes from rare, ex- 
tended regions which can host gapless excitations, while the 
typical behavior (captured at the mean-field level) is that of a 
gapped system. 



B. Two-species Bose glass 

As seen in the previous subsection, the DLM phase and the 
BG phase are indeed continuously connected, and the DLM 
phase has a nature which is analogous to that of the BG phase. 
In Section lTV Fl we pointed out that the DLM phase is charac- 
terized by exponentially localized quasiparticles correspond- 
ing to LMs antiparallel to the field, while the BG phase is char- 
acterized by the exponential localization of dimer triplets. For 



0.275 < J'/J < 0.45 localized quasiparticles of hoth species 
appear above the transition curve for the clean system (which 
is the condition for the presence of dimer triplets). The ap- 
pearance of localized dimer triplets before the LM quasipar- 
ticles localized on LM clusters disappear upon increasing the 
field guarantees the gapless nature of the many-body spectrum 
throughout the quantum disordered regime. Hence the phase 
diagram region for J'/ J > 0.275 comprised between the two 
curves /ibg( J'/-^) for 'he doped system and hc{J' / J) for the 
clean system represents a two-species Bose glass (2BG). 

Fig. [13] shows the magnetization curve along a section of 
the phase diagram (J'/J = 0.35) crossing the 2BG phase. 
Similarly to the case of weaker inter-dimer coupling thor- 
oughly investigated in this paper, clear PPs, corresponding to 
LM dimers and trimers, are observed at rriu/mu^^^ ~ 0.73 
and rriu/mu^^^ ~ 0.95, where rriu^^^ = p/2 is the satura- 
tion magnetization of the LMs. The plateau at m„ ' is on 
the contrary completely removed by the appearance of dimer 
triplets at a field hc{J' /J) ~ 0.4J, which marks the onset of 
the 2BG phase. 




FIG. 13: Field evolution of the uniform magnetization mu and 
the uniform susceptibility Xu in the site-diluted dimer systems with 
J' /J = 0.35 and p — 1/8. Although there is no true magnetization 
plateau, two pseudo-plateaus rriu/m'^^ ~ 0.73 and mu/m^^ ~ 
0.95 are resolved. 



C. Quantum Scaling and Critical Exponents 

Next we investigate the properties of the quantum critical 
line separating the quantum disordered phase from the or- 
dered phase in the phase diagram of Fig. [12] As already 
pointed out above, for weak inter-dimer coupling (J' ^ J) 
this critical line divides ordered and disordered phases which 
have a different nature when considering the low-field or 
the intermediate-field region. At low fields, the critical field 
(/iDLivi) divides the OBD phase from a DLM phase, while at 



13 



higher fields, the critical field (/ibg) divides a BG phase of 
dimer triplets from a condensate of those triplets (correspond- 
ing to an XY ordered antiferromagnetic phase). Hence a fun- 
damental question arises whether the quantum critical curve 
is characterized by universal critical exponents or not. 
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FIG. 14: Finite size scaling of the correlation length in the vicinities 
of the critical fields ^dlm (left) and /ibg (right). The inter-dimer 
coupling is J' — J/4, and the dilution concentration is p = 0.125. 



We answer this question by a detailed numerical scal- 
ing analysis of the transition line. The correlation length is 
extracted from the transverse structure factor Eq. Q by a 
second-moment method,-'^-^" assuming that S-^{k) follows a 
Lorentzian shape close to fc = (tt, tt): 




S'^(7r,7r) 



(15) 



where Afc 



2tt 
L 



S^iTT + Ak,7r) 
The scaling form of the correlation length is 

^ = Lf^[L^/''{h-h,)], (16) 



where f^[-] denotes the corresponding scaling function. This 
defines the correlation exponent ly. The scaling plots of the 
correlation length for J'/ J = 1/4 around the two critical 
fields huLM = 0.007(1) and /ibg — 0.67(1) are shown in 
Fig- [HI and they both show the best collapse for an exponent 
i/« 1.0(1). 

To further characterize the universality class of the transi- 
tion we consider the scaling properties of the staggered mag- 
netization and spin stiffness, whose scaling forms are 

m, = L-^/VrnJi'/"^(/i-M], (17) 
Ps = L-^''^^-'^fpAL'/''ih-K)]. (18) 

Here /3 is the order parameter exponent, z is the dynamical 
exponent, and 13 = 2 is the spatial dimension of the system. 
Fig. [15] shows the associated scaling plots, which give best 
collapse for critical fields /idlm and /ibg and for a correla- 
tion exponent all consistent with the values estimated via 
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FIG. 15: Finite size scaling of the staggered magnetization and the 
spin stiffness in the vicinities of the critical fields /idlm (left) and 
/iBG (right). The model parameters are the same as in the caption of 
Fig.ini 



the scaling of the correlation length. Moreover we extract the 
other exponents (3 = 0.9(1) and z = 2.0(1) at both /idlm and 

hBG- 

For a BG-SF transitions, it was theoretically predicted that 
z = D = 2*^ in agreement with our results. Moreover for a 
well defined transition to occur in a disordered system the cor- 
relation length exponent ly' must satisfy the Harris criterion^ t 
ly ^ 2/D = 1, which is also verified by our results. Repeating 
the same scaling analysis along the critical line of the phase 
diagram we obtain analogous estimates for the critical expo- 
nents, as shown explicitly in Fig.[T6]for the case J'/ J ~ 0.35; 
this is indeed a non-trivial result, as the BG phase close to /ibg 
in this case is a 2BG, as discussed in Section IVBl and across 
the transition at /ibg the dimer triplets condense, while the 
LM quasiparticles remain localized. 

Hence all the above results point towards a critical curve 
with universal critical exponents defining the SF-BG quantum 
phase transition in 13 = 2. The exponents [3, v and z are fully 
consistent with those estimated in Ref. 32 for the SF-BG tran- 
sition in a differently coupled-dimer system (a strongly cou- 
pled bilayer) with site dilution, and therefore we argue that 
they apply to the condensation transition of a generic two- 
dimensional dirty-boson system with incommensurate parti- 
cle density. 
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FIG. 16: From top to bottom: finite size scaling of correlation length 
^, staggered magnetization rus, and spin stiffness ps at critical fields 
ftDLM = 0.05(1) (left) and /ibg = 0.56(1) (right) in the site-diluted 
coupled dimers with J'/ J = 0.35 andp =1/8. 



VI. CONCLUSIONS 

Making use of Stochastic Series Expansion Quantum 
Monte Carlo, we systematically investigated the zero- 
temperature phase diagram of weakly coupled dimers of S* = 
1 /2 spins with site dilution and in the presence of a magnetic 
field. We focused our attention on the specific case of a planar 
array of coupled dimers. In particular we find that the phase 
diagram is much richer than the clean system, and also richer 
than what has been predicted at the mean-field level^^. In par- 
ticular, we show that the impurity-induced ordered state of lo- 
calizes moments (LMs) at zero field is destroyed by a tiny field 
for sufficiently small dilution, and that the field drives the sys- 
tem to a novel, disordered-local-moment (DLM) phase, which 
extends over a very large field range. We elucidate the micro- 
scopic nature of this gapless and quantum-disordered phase 
by comparing results for the diluted coupled-dimer array with 
those obtained for a network of S* = 1/2 LMs interacting 
via effective couplings exponentially decaying with the dis- 
tance. We observe that the DLM phase is characterized by the 
fact that spins pointing oppositely to the field remain local- 
ized on rare clusters of nearest-neighboring LMs. A spin-to- 
boson mapping shows that the DLM phase is akin to a Bose 
glass phase of localized hardcore bosons. In the DLM phase. 



the magnetization curve of the system shows a hierarchy of 
pseudo-plateaus marking each the full polarization of a class 
of LM clusters (dimers, trimers, etc.). The specific quantum 
nature of the DLM phase is revealed by direct comparison 
with data for the classical limit of the LM network. 

A systematic study of the phase diagram of the site-diluted 
planar array of dimers, as a function of the ratio of inter- 
dimer coupling to intra-dimer coupling and of the strength 
of the applied field, reveals that the DLM phase is continu- 
ously connected to the Bose glass phase occurring at higher 
fields, and associated with the Anderson localization of triplet 
quasi-particles appearing on the intact dimers. This continu- 
ous disordered and gapless phase separates the gapped phase 
of fully polarized LMs from the ordered phases - either the 
impurity-induced phase present also at zero field, or the field- 
induced ordered phases characterized by condensation of the 
triplet quasi-particles. This special feature of the phase dia- 
gram is consistent with what expected for dirty -boson systems 
on a lattice. Moreover the critical exponents for the quantum 
phase transition separating the gapless disordered phases and 
the ordered phases appear to be same regardless the specific 
nature of the ordered and disordered phases which are con- 
nected by it. This signals the existence of a novel universal- 
ity class for the Bose glass-to-superfluid transition in 2D with 
V w 1.0, z w £) = 2, and /3 w 0.9. 

The above findings show that the observation of unconven- 
tional disordered phases without a gap, induced by disorder, is 
quite realistic in the context of doped spin-gap antiferromag- 
nets. In fact, recent studies have focused on the properties 
of the bond-disordered spin ladder IPA-Cu(Clj, Bri_j,)3^^"^^, 
giving the first evidence of a magnetic Bose glass in this 
system^^*^. In the context of bond disorder the Bose glass 
appears as a consequence of localization of triplets induced 
by intense fields^^, while no DLM phase is expected. On 
the other hand, site dilution of the magnetic lattice, realized 
by non-magnetic doping, would realize a much broader Bose 
glass region (comprising the DLM phase), to be observed also 
at small magnetic fields. This is not at all a negligible aspect: 
indeed the smallness of the field for which the DLM phase sets 
in (which can be two or more orders of magnitude smaller than 
the dominant antiferromagnetic coupling) makes the DLM 
phase observable also in compounds which have a large spin- 
gap (of the order of 10 K or higher), and for which the physics 
of spin-triplet condensation is not easily accessible in exper- 
iments. A possible candidate compound is the recently in- 
vestigated spin ladder compound Bi(Cui_:rZnj;)2P06r^ fea- 
turing a spin gap of ^ 35 K with zero doping {x = 0). 
Other potential candidates include the spin-ladder compounds 
(C5Hi2N)2CuBr4^'^ '^ and IPA-CuCla^^ with Zn of Mg dop- 
ing of the Cu ions. 
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